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We report on the experimental demonstration of stimulated Raman adiabatic passage (STIRAP) 
in a Tm 3+ : YAG crystal. Tm 3+ :YAG is a promising material for use in quantum information process- 
ing applications, but as yet there are few experimental investigations of coherent Raman processes 
in this material. We investigate the effect of inhomogeneous broadening and Rabi frequency on 
the transfer efficiency and the width of the two-photon spectrum. Simulations of the complete 
Tm 3+ : YAG system are presented along with the corresponding experimental results. 
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I. INTRODUCTION 

Rare earth ions doped into inorganic crystals have re- 
cently been considered for use in the field of quantum 
information processing [l], 0, H, 0, H, ■ These materials 
are promising as they offer very long coherence times in 
both the hyperfine j3] and optical transitions 0. Cur- 
rently the majority of quantum information processing 
demonstrations have been performed in crystals doped 
with either praseodymium or europium. Unfortunately 
these materials can only be optically addressed with dyc- 
lasers. The task to stabilize these lasers to the sub-kHz 
lincwidths required for quantum information processing 
experiments is so challenging that only a few groups in 
the world have achieved such laser sources 0, [l(| [ll[ . 

Another disadvantage of praseodymium and europium 
doped crystals is the small hyperfine splitting, which is 
of the order of ten's of MHz. This small separation lim- 
its possible quantum information storage demonstrations 
in these crystals as the bandwidth of existing quantum 
sources does not match this narrow hyperfine splitting. 
It is therefore worthwile to investigate other rare earth 
ions which do not suffer these disadvantages. 

One such promising rare earth ion is thulium. The 
wavelength of thulium is 793 nm which falls in the range 
of easily stabilized semi-conductor lasers. The only iso- 
tope of thulium has an / = 1/2 nuclear spin but un- 
fortunately does not exhibit hyperfine structure due to 
J quenching [r|. It has proven possible to create a 3- 
level A-system in this material via the application of a 
magnetic field with a particular orientation [l3[ . The ad- 
vantage of this A-system is that it possesses an adjustable 
ground state hyperfine splitting via the magnitude of the 
applied magnetic field. Thus Tm 3+ :YAG does not suffer 
from a limited bandwidth for possible quantum informa- 
tion applications and is optically accessible with easily 
stabilized lasers. 

There has been a startling lack of investigation into Ra- 
man coherent processes in this material. The first nuclear 
spin coherence via optical excitation has only recently 
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been demonstrated [14j and it is certainly obvious that 
these types of processes need to be studied further if this 
material is to be useful for quantum information process- 
ing applications. Stimulated Raman adiabatic passage 
(STIRAP) [UGl 

is one such coherent process which 
deserves attention in this material. It provides complete 
coherent population transfer in a A-system without suf- 
fering from radiative losses and is robust against slight 
pulse variations. STIRAP is thus important in the field 
of quantum information processing for use in the prepa- 
ration and manipulation of ions. 

In this paper we will first discuss the theory of STI- 
RAP using the simple 3-level A-system at two-photon 
resonance. We will then present numerical simulations 
of the full Tm: 3+ :YAG level structure alongside simula- 
tions of the simple 3-level A-system. Following this will 
be a presentation of experimental and numerical results 
along with discussion and conclusions. 



II. THEORY 

For a simplified discussion of STIRAP we consider the 
3-level A-system shown in Fig. [TJ The objective is to 
efficiently transfer population from level |1) to level |3) 
via the use of two temporally overlapping laser pulses. 

The ground state hyperfine levels are split by A ff . The 
detuning between the Stokes pulse and optical resonance 
is denoted A and is referred to as the optical or onc- 
photon detuning. The detuning of the pump pulse from 
two-photon resonance is denoted 8 and corresponds to the 
two-photon detuning. In the simplified case we assume 
that 6 = 0, ie. we are at two-photon resonance. The 
system is initially prepared with all the population in 
level |1) and two, non-simultaneous, laser pulses are ap- 
plied. The Stokes pulse, applied first, excites the, initially 
empty, transition |3) — |2) and the pump pulse, appplicd 
temporally overlapping the Stokes pulse, excites the tran- 
sition |1) — 1 2). The pump and Stokes pulses have the fol- 
lowing, respective, Rabi frequencies: S! p = ~^i2£ P (t)/h 
and fig = —H32£s(t)/H. The dipole transition moments 
are the /Uy's and the applied, time-dependent, electric 
fields of the laser pulses are the £ SjP (i)'s. 

The eigenstates of the instantaneous RWA Hamilto- 
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FIG. 1: Three level A-system where the ground state hyper- 
fine levels are split by A a . The optical, or one-photon, de- 
tuning is labelled A and is the frequency difference between 
the Stokes pulse and the transition |3)-|2). The two-photon 
detuning is labelled S and is the frequency difference between 
the pump pulse and two-photon resonance. The Rabi fre- 
quencies of the pump and Stokes pulses are labelled Sl p and 
fi s respectively. 



nian of this system are well known and are given by the 
following linear combinations of the bare states |1), |2) 
and 1 3): 

\a + ) = sin sin ^> |1) + cos</> |2) + cos0sin</> |3) 
\a°) = cos0|l) -sin0|3) 

|a~) = sin0cos</) |1) — sin0 |2) + cos0cos(/> |3) (1) 

where the time-dependent mixing angles, 9 and <f>, are 
given by: 



tan 9 



tan 2(f) 



(2) 



and we have the following time-dependent, dressed state 
eigenvalues: 

lj° = 



ur^t) = A/2±l/2^A 2 + n 2 p (t) + fl 2 8 (t) (3) 

After examining the eigenstates it is obvious that two 
dressed states, l^), contain the 'leaky' bare state |2). 
Thus these dressed states cannot be used for complete 
population transfcrral as some population will be pro- 
moted to the excited level, |2), where it will experience 
radiative losses and therefore be lost to the final state 
13). 

Inspection of the dressed state \a°), known as the 
'dark' state as it contains no component of the electronic 
excited state, demonstrates that it should be possible to 
perform efficient population transfcrral using this state. 
If the Stokes pulse is applied before the pump pulse then 



at time t —> — oo we have the mixing angle = 0° and 
the dressed state |a°) corresponds to the bare, populated, 
state |1). At the end of the interaction, at time t — ► oo, 
where the pump pulse has been applied after the Stokes 
pulse, we have a mixing angle of 9 = 90° and the state 
| a ) corresponds to the bare state |3). 

This implies that provided the evolution of the sys- 
tem is adiabatic and that the Stokes pulse preceedes the 
pump pulse the system will remain in the |a°) state, or 
the 'dark' state, throughout the interaction. Hence the 
population initially in level |1) will be transferred into 
level 1 3) whilst never populating level |2). Thus the pop- 
ulation transfcrral does not depend on either the optical 
coherence lifetime or the decay rate from the upper level 
provided the adiabacity condition is satisfied. 

The adiabaticity condition requires the energy spacing 
between the eigenvalues to be much larger than the dy- 
namic coupling term, given by 0(t), where the overdot 
represents the time derivative. The adiabacity condition 
is US: 



± 



v /A2 + ^(t) + ^ w 



(4) 



Thus for any given shape of the laser pulses, £l p and £! s , 
the adiabacity condition can be evaluated at any time, t. 

At optical resonance this adiabacity condition is gen- 
erally satisfied, for smooth pulses, if: 



SleffT > 1 



(5) 



where t is the pulse duration and the effective Rabi fre- 
quency is: 



eff 



(6) 



Far from optical resonance the adiabatic condition, for 
smooth pulses, becomes instead: 



ni ffT » 



Al 



(7) 



which sets both an upper limit on the optical detun- 
ing and a lower limit on the Rabi frequencies for which 
population transfcrral is efficient. Therefore carefully se- 
lected pump and Stokes pulses, tuned to two-photon reso- 
nance, will transfer population adiabatically over a spec- 
tral range given by the adiabacity condition of Eqn. [7] 

The eigenvectors and eigenvalues of the Hamiltonian 
involving a non-zero two-photon detuning, 5, are much 
more complicated than those of the two-photon reso- 
nance case. These eigenvectors and eigenvalues are ex- 
plicitly solved in Ref. [l8| in which it is demonstrated 
that all three eigenvectors contain a contribution of the 
excited state level. It was also shown that away from 
two-photon resonance there is no adiabatic state which 
connects the initial bare state at the beginning of the se- 
quence with the desired final bare state at the end of the 
interaction [lc| . This implies that adiabatic processes are 
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not the means for population transfcrral in this type of 
system. When two-photon detuning is present the major- 
ity of the population is instead transferred via diabatic 
processes which involve transisitions between the adia- 
batic states. 



III. SIMULATIONS 

Simulations of the Tm 3+ :YAG system were performed 
in order to aid our understanding of the system. The full 
Tm 3+ :YAG level scheme is actually quite simple when 
compared to that of the full praseodymium and europium 
systems. Praseodymium and europium both have a nu- 
clear spin of / = 5/2 and thus when in a magnetic field 
the level scheme consists of six energy levels in both the 
ground and excited electronic states. Thulium on the 
other hand has a nuclear spin of / = 1/2 and therefore 
when it is placed in a magnetic field the energy level 
structure is much simpler, consisting of only two Zceman 
energy levels in both the ground and excited electronic 
states. 

This simple energy level structure of Tm 3+ :YAG, dis- 
played in Fig. [21 is thus easy to completely model, even 
when including the additional metastable energy level. 
In the full Tm 3+ YAG level structure the ground and 
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FIG. 2: Level structure of Tm 3+ :YAG in a magnetic field, the 
ground and excited electronic states are split into two Zee- 
man levels which are labelled with their nominal nuclear spin 
numbers. The ground and excited state hyperfine splittings, 
labelled A 9 and A e respectively, are dependent on the mag- 
nitude of the applied magnetic field, with A e = A g /2.5. The 
two-photon detuning and the Rabi frequencies of the tran- 
sitions are as described in the 3-level system. The optical 
detuning is the frequency difference between the Stokes pulse 
and the transition |3)-|4). 

excited state splittings arc labelled as A g and A e respec- 
tively, and are dependent on the magnitude of the applied 
magnetic field, with A e = A g /2.5. The optical detuning, 
A, and the two-photon detuning, <5, are as defined in the 
previous section. 

One complication to the Tm 3+ YAG model arises from 
the fact that the Stokes and pump pulses can excite all 



four optical transitions. This is possible due to the lack 
of polarisation selectivity of the optical transitions and 
the mixing of the nuclear spin states, which relaxes the 
spin selection rule, Am = 0. This relaxation results in 
the spin-flip transitions, |3)-|2) and |1)-|4), having a tran- 
sition dipole moment of 0.37 times that of the transitions 
not involving a spin- flip [l9j . The nominal nuclear spin 
labels in Fig. [2] thus represent the strong and weak tran- 
sitions. 

Whilst the Stokes and pump pulses are applied with 
equal intensities the Rabi frequencies of the transitions 
are different depending on which transition is being ex- 
cited. The Rabi frequency of the weak transition is 0.37 
times the value of the Rabi frequency of the strong tran- 
sition. In all the following discussions the Rabi frequency 
listed is that of the strong transition. 

The metastable state in Fig. [5] is included as a kind 
of reservoir in the model. This state has a much longer 
radiative lifetime, of 10 ms [2(|, compared to that of the 
optical excited level, of 800 /is [2l|. In addition to this, 
most of the population in the optical excited level decays 
to the metastable state rather than the ground state, in 
a ratio of 3:1 [13]. 

The optical Bloch equations for this full energy level 
scheme, including the metastable state and the fact that 
multiple transitions are excited by the pump and Stokes 
pulses, are numerically solved in the interaction picture. 
We have also included population and coherence decay 
terms in the equations and the transfer efficiency is de- 
termined 300 /xs after the application of the pump pulse. 
The Stokes and pump pulses both have a Gaussian line- 
shape with a FWHM (of intensity) of 30/-\/2 /is and a 
pulse delay of -17 /is, with the Stokes pulse preceding 
the pump pulse. 

In Fig. [3] we compare the optical detuning scans ob- 
tained from the 4-level model described above and the 
standard 3-level model of the level structure shown in 
Fig. [TJ This 3-lcvcl model is simply the numerical solu- 
tion to the standard optical Bloch equations where we 
have included dephasing and relaxation terms and the 
fact that we have a strong and weak transition. Both 
these simulations were performed at two-photon reso- 
nance, ie. at S = 0. In this case the transfer efficiency is 
determined from the population difference between the 
excited level, level |4) in the 4-level model, and level |3). 

Looking at only the one-photon linewidth centred at 
A = 0MHz in Fig. [3] for the 3-levcl and 4-level models 
we see that the introduction of the fourth level and the 
multiple excitations of the transitions results in an asym- 
metry to the one-photon lincshape as well as reducing the 
overall linewidth. The one-photon linewidth (FWHM) of 
the 4-level model is 1.3 MHz as compared to 1.5 MHz for 
the 3-level model. From the particular pulse parameters 
used here the adiabacity criterion of Eqn. [7] reduces to 
6.3 MHz ^> |A|. The linewidths of both the models agree 
with this criterion but this equation does not explain the 
asymmetry to the 4-level model. 
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FIG. 3: (colour online) The optical detuning scan from the 4- 
level model, dashed line (red), and the standard 3-level model, 
solid line (black). In this case the excited state splitting was 
2.84 MHz, the Rabi frequency on the strong transition was 
510 kHz, the pulse duration was 30/^/2 (is and the pulse delay 
was -17 /is. These models were both performed at two-photon 
resonance and we record the transfer efficiency at each optical 
detuning value, A. 

In Fig. [3] we also see the contributions from both A- 
systems in the 4-level model. At A = MHz we are 
on-resonance with the excited state level |4) whilst at 
A = 2.48 MHz we are now on-resonance with the other 
A-system, with the excited state level |2). It can be seen 
in this figure that the maximum transfer efficiency is the 
same for both A-systems even though the Rabi frequen- 
cies of the Stokes and pump pulses have been switched. 
It is thought that the transfer efficiency dramatically de- 
creases at A = 1.24 MHz, the midway point between the 
two excited state levels, due to an interference effect be- 
tween these two A-systcms. It is this interference effect 
which causes the asymmetry and reduced linewidth of 
the 4-level model over the 3-level. 

In Fig. [4] is shown the transfer efficiency determined 
from the 4-level model as both the one- and two-photon 
detunings are varied. The effect of the second excited 
state hyperfine level can be seen in this figure. The ef- 
ficient population transfer centered at A = is due to 
the pump and Stokes pulses being optically resonant with 
level |4) of Fig. [5] Whilst the efficient transfer centered 
at A = 2.84 MHz is due to the pump and Stokes pulses 
being applied nearly resonant with level |2) of Fig. [21 the 
second excited state hyperfine level. The lack of popu- 
lation transfer for optical detuning values between these 
two levels is thought to be due to the interference effect 
of the two A-systems. 

We also compare the standard 3-level model and the 
4-level model described above for a two-photon detuning 
scan. In this case we perform the simulations at optical 
resonance, ie. at A ~ 0, whilst we vary the frequency 
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FIG. 4: The transfer efficiency determined from the 4-level 
model for different values of the two-photon and optical de- 
tunings for an excited state splitting of 2.84 MHz, a Rabi 
frequency of 300 kHz and a pulse delay of -17 /xs. 

difference between the pump pulse and two-photon reso- 
nance, 5. 
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FIG. 5: (colour online) Two-photon detuning scan at optical 
resonance, A = 0. The solid line (black) is the simulation 
of the standard 3-level model and the dashed line (red) is 
the result from the 4-level model. In these models the Rabi 
frequency on the strong transition was 400 kHz, the ground 
state splitting was 7.1 MHz, the excited state splitting was 
2.84 MHz and the pulse delay was -17 [is. 

We see in Fig. [5] that the addition of the fourth level 
makes less difference to the two-photon detuning scan 
than it did for the optical detuning scan. The two- 
photon linewidth is again slightly narrower for the 4-level 
model. The two-photon linewidth (FWHM) of the 3-level 
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model is 360 kHz whilst that of the 4-level is 320 kHz 
and there is again a slight asymmetry to the 4-level line- 
shape. Adiabatic population transferral can only occur 
at two-photon resonance. Thus, we see in Fig. [5] that 
there is a large amount of population transferral due to 
non-adiabatic, or diabatic, transitions in both models. 
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IV. EXPERIMENTAL SET-UP 

When performing STIRAP experimentally the system 
is excited with an extended cavity semi-conductor laser 
operating at 793 nm. It is stabilized with a high-fincssc 
Fabry-Perot cavity via a Pound-Drever-Hall scrvoloop 
with a frequency stabilisation of 200 Hz over 10 ms [23| . 
The laser is then amplified with a semi-conductor tapered 
amplifier (Toptica BoosTA). The experimental set-up is 
shown in FigO AOMs 1 and 2 are arranged in a double- 
pass configuration so that the light on the sample is not 
spatially shifted with a frequency shift. These AOMs arc 
driven by a dual-channel 1 Gigasample/s arbitrary wave- 
form generator (Tektronix AWG520) which can produce 
arbitrary amplitude and phase shaping. 

After passing back through the polarizing beam- 
splitter a common polarization is given to both compo- 
nents via a Glan prism. This recombined beam is then 
coupled into a 2 m long single-mode fibre. The light po- 
larization direction is optimized with a half-wave plate 
to ensure maximum Rabi frequency. This beam is then 
focussed onto the 0.1 at.% Tm 3+ :YAG sample which is 
5 mm in the light propagation direction. The sample is 
cooled to 1.8 K in an Oxford 6T Spectromag SM4 cryo- 
stat. 

After the cryostat the beam passes through a 50 fj,m 
pinhole and is then imaged onto an avalanche photodiode 
(Hamamatsu C5460). The Gaussian shaped cross-section 
of the laser beam means that those ions positioned at the 
edge of the beam diameter do not experience the maxi- 
mum Rabi frequency. The pinhole ensures only those ions 
in the centre of the beam and thus subjected to the max- 
imum Rabi frequency of the pump and Stokes pulses are 
probed. A final AOM protects the detector from strong 
light pulses. 

The applied magnetic field, generated by supercon- 
ducting coils, is applied in the direction optimizing the 
branching ratio [13| . The magnitude of this applied mag- 
netic field determines the splitting of the hyperfine levels. 
The ground state splitting is given b y 41 MHz/T with an 
excited state splitting of 16 MHz/T [lj]. The magnitude 
of this applied field was of the order of 0.2-0.5 T. 

The inhomogencous width of the optical transition is 
^25 GHz which is much larger than the hyperfine level 
splitting. Thus a single laser frequency can drive transi- 
tions in different ensembles of ions at the same time. In 
order to drive and probe the STIRAP process in a single 
ensemble of thulium ions with a well-defined transition 
frequency we perform a preparation sequence based on 
optical hole burning. 



FIG. 6: Experimental set-up of the system, acousto-optic 
modulators are labelled as AOM. AOMl and 2 are arranged 
in a double-pass arrangement so the light on the crystal is not 
spatially shifted when the frequency changes. These double- 
pass AOMs are centred on 110 MHz (AA OptoElectronics) 
and provide the pump, Stokes, probe and preparation pulses. 
AOM3 is used as a gate to protect the detector from the in- 
tense preparatory and Gaussian-shaped pulses. 



The preparation sequence consists of repeated burning 
over a range of 3 MHz via scanning the laser. This optical 
pumping creates a trench in the absorption profile. After 
this preparation sequence the pump and Stokes pulses are 
applied. The frequency of the Stokes pulse is chosen such 
as to place it in the centre of this created transmission 
window and thus initially it experiences no absorption. 
All the population is initially in level |1) of Fig [21 which 
is a necessary requirement for STIRAP. 

The pump pulse is applied at a frequency outside this 
transmission window and the frequency difference be- 
tween the Stokes and pump pulses matches the ground 
state splitting of the system, A 9 . After STIRAP has 
modified the population distribution the absorption of a 
weak probe pulse provides a measurement of the popula- 
tion difference between the final level |3) and the optical 
excited state. 

The temporal profiles of the Stokes and pump pulses 
are nominally Gaussians, but there was a saturation ef- 
fect from one of the AOMs with altered the shape of the 
Stokes pulse thus changing the Gaussian profile. The 
FWHM of intensity for the two pulses was nominally 
30/^/2 /is. The delay between the pump and Stokes 
pulses is variable, but is -17 fis, unless otherwise stated, 
with the Stokes pulse preceeding the pump pulse. The 
probe pulse has a rectangular profile with a pulse dura- 
tion of 1 ms, during which the frequency is scanned over 
20 MHz, and is applied 300 /is after the pump pulse. In 
the interaction region of the crystal the pump, prepa- 
ration and Stokes beams all have a beam diameter of 
110 /im. 

Shown in Fig. is a typical experimental spectrum 
recording the transmission of the weak probe pulse as its 
frequency is scanned. The overlapping red line is the ex- 
pected probe transmission for this particular optical den- 
sity after the preparation and STIRAP pulse sequences 
have been applied. In this case the ground state splitting 
was A g = 4.55 MHz and the excited state splitting was 
A e = 1.82 MHz. 

The various trenches of either increased or decreased 
transmission are related to the initial preparation step 
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FIG. 7: (colour online) A typical experimental scan of the 
transferred population where the probe frequency, v s is 

the Stokes frequency, v v is the pump frequency and A e is the 
excited state hyperfine splitting. The ground state splitting 
was A g = 4.55 MHz with an associated excited state splitting 
of A e = 1.82 MHz. The Rabi frequency of the strong tran- 
sition was 1 MHz. The black line is the experimental output 
whilst the red line is a fit to the experimental scan for this 
particular optical density. 

of the sequence. After the preparation sequence there 
is increased transmission around the Stokes frequency, 
zAs, associated with the population having been optically 
pumped into the other ground state hyperfine level. Cor- 
respondingly there is decreased transmission around the 
pump frequency, v p , associated with the increased popu- 
lation in this level. The decreased transmission also seen 
around v s — A g is due to the symmetrical nature of the 
holcburning spectrum. 

The preparation sequence actually creates two ensem- 
bles of thulium ions participating in the STIRAP se- 
quence. One ensemble is that displayed in Fig. [5] where 
the Stokes pulse is applied on the strong transition. The 
other ensemble is the opposite case, where the Stokes 
pulse is nearly resonant with the weak transition. In 
each ensemble the situation is the opposite for the pump 
pulse. 

The STIRAP sequence of the temporally overlapping 
Stokes and pump pulses modifies this initial population 
distribution and the associated STIRAP absorption and 
transmission features are observed overlaid on the prepa- 
ration trenches. 

The main STIRAP absorption feature is observed at 
the Stokes frequency which contains contributions from 
both ensembles of ions. A STIRAP absorption feature is 
also seen at the Stokes frequency plus the excited state 
splitting, v s + A e , which contains contributions from the 
ensemble where the Stokes pulse interacts with a strong 



transition. There should also be a STIRAP absorption 
feature at v s — A e , but the contributions to this feature 
only occur from the ensemble where the Stokes pulse in- 
teracts with a weak transition and thus is much smaller 
than the features at v s and v s + A e . 

The corresponding STIRAP transmission features are 
observed at the pump frequency and at the pump fre- 
quency minus the excited state splitting, v p — A e . Again 
the absorption feature which should be observed at 
v p + A e contains contributions solely from the ensemble 
where the pump pulse interacts with a weak transition. 

The important intensity levels are labelled Ii, I2 and 
I3 in Fig. Ii corresponds to equal population in the 
two hyperfine ground state levels, |1) and |3), I2 corre- 
sponds to the amount of transferred population in level 
1 3) after the STIRAP sequence, whilst I 3 corresponds to 
zero population in the final level |3). With these three 
intensity levels it is thus possible to determine the trans- 
fer efficiency using Beer's law. In this case the transfer 
efficiency was ^94%. 

The level scheme in Fig. [5] demonstrates how our 
method of optical detection functions. The STIRAP 
experiment performed here is actually a two step pro- 
cess. First we perform STIRAP. The two-photon de- 
tuning and the ground state splitting are related to the 
pump-Stokes detuning by A g — 8 = v v — v s . The ground 
state splitting is fixed by the applied magnetic field am- 
plitude. However, due to the inhomogeneous width of 
the ground state hyperfine transition in our system, the 
external field can only set an average value, (A g ). There- 
fore, the given pump and probe frequencies can only spec- 
ify an average value of the two-photon detuning, given by 
(6) = (A g ) + v s - v p . 

In the second step we scan the frequency of the weak 
probe pulse, v pr . The probe field precisely selects the ions 
with a definite value of the STIRAP optical detuning, A, 
given by A = Vp r — v s . Indeed the probing spectral resolu- 
tion is given by the homogeneous width, which amounts 
to a few kHz. When v pr = v s we are at optical resonance 
and the maximum transfer efficiency is observed. 

In Pr 3+ :Y2Si05 the inhomogeneous broadening of the 
hyperfine transition is small with respect to the prac- 
tical Rabi frequencies. Therefore, as far as STIRAP is 
concerned, the two-photon detuning is precisely defined 
by the pump and Stokes frequency difference, v v — v s . 
Specifically, appropriate adjustment of this frequ ency 
difference can lead to two-photon resonance [2J, |25I |. 
In Tm 3+ :YAG, the Raman inhomogeneous width grows 
with the magnitude of the applied magnetic field at the 
large rate of ^400 kHz/T. Hence the incidence of this in- 
homogeneous broadening on the STIRAP efficiency has 
to be examined. 



V. RESULTS AND DISCUSSION 

We will first present initial investigations of STIRAP in 
Tm 3+ :YAG, characterising the transfer efficiency and the 
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width 

FIG. 8: Experimental level scheme, this is the simplified 
scheme of that shown in Fig. [5] where we have removed level 
1 2} for ease of reference. The Stokes and pump pulses are 
optically resonant with a particular A scheme comprising of 
levels |1), 1 3) and The probe pulse is then scanned and 
interacts with different A-systems with optical detunings of 
A, determined during the STIRAP step. The ground state 
splitting and two-photon detuning values, (A 9 ) and (5), are 
average values due to the Raman inhomogeneous braodening 
present in our system. 
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FIG. 9: (colour online) We vary the delay between the Stokes 
and pump pulses. The dashed line (red) is the delay scan 
of the standard 4-level model whilst the dash-dot line (blue) 
is the 4-level model where we have included the inhomoge- 
neous width to the Raman transition, henceforth to be la- 
belled 4-level+, and the squares are the experimental data. 
The ground state splitting was A 9 = 4.55 MHz, the inho- 
mogeneous width of the Raman transition was ~50 kHz, the 
Rabi frequency of the strong transition was 1 MHz and the 
pulse widths were 30/\/2 fis. The simulations were performed 
at both optical and two-photon resonance. 



one- and two-photon detuning scans. Following these ini- 
tial STIRAP investigations we will present numerical and 
experimental data demonstrating the effect of the inho- 
mogeneous width and Rabi frequency on both the maxi- 
mum transfer efficiency and the two-photon lincwidth of 
this system. 

In our first experiment we varied the delay between 
the Stokes and pump pulses, covering both the Stokes 
preceding the pump, negative delays, and the pump pro- 
ceeding the Stokes, positive delays. The resulting trans- 
fer efficiencies were monitored for these different pulse 
orderings and are the squares displayed in Fig [5] 

We see that for optimum negative delays the experi- 
mental transfer efficiency reaches nearly 100%. Signifi- 
cant population transferral is also observed for positive 
delays, where the pump precedes the Stokes pulse. In 
both these cases the population transferral is significantly 
larger than that transferred when the pump and Stokes 
pulses arc applied coincident. The maxima in the ex- 
perimentally measured transferred population show the 
characteristic plateaus which are a signature of the STI- 
RAP process. These plateaus imply that the transferral 
process is robust against small pulse variations. 

In Fig. [5] we see that the maximum transfer efficiency 
occurs when the pulse separation of the Stokes and pump 
pulses is equal to the pulse width, at delay~-21 fis. It 
has been previously discussed [HI, [26| that the maximum 
transfer efficiency is to be expected at this pulse separa- 



tion for Gaussian shaped pulses as the adiabacity condi- 
tion of Eqn. 21 taking A = for simplicity, is best satis- 
fied throughout the interaction at this pulse delay. Thus 
we have the lowest losses via non-adiabatic coupling to 
the 'leaky' dressed states la 1 * 1 ) and the maximum transfer 
efficiency. 

The significant population transferral for pos itive de- 
lays is consistent with 6-STIRAP [H H3, III]. In b- 
STIRAP the system remains in the bright state and the 
excited state is populated during the interaction. Radia- 
tive losses are therefore possible resulting in a reduction 
in the transfer efficiency. 6-STIRAP is only observed in 
systems where the excited state radiative lifetime is long 
compared to the interaction time. In Tm 3+ :YAG the 



lifetime of the excited state is 800 fjs [21( , which is much 
longer than the interaction time of the pulses. 

We see that the experimental data in Fig. [9] agrees 
well with both the standard 4-level and the 4-level model 
where we have included the inhomogeneous broadening 
to the hyperfine transition, henceforth to be known as 
the 4-level+ model. In performing these simulations and 
in subsequent modeling we have included the saturation 
effect of the AOM in our models. 

When including the significant inhomogeneous 
linewidth of the hyperfine transition in our 4-level+ 
modcl we assume the inhomogeneous line is a Gaussian 
centred on (8). Our position in the inhomogeneous line 
is now specified as 5 and thus for each value of (8) we 
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must integrate over these S's, whilst also including the 
Gaussian lineshapc. The two-photon detuning profile 
in Fig. [5j for the standard 4- level model, is the type of 
profile we use to perform this integration. In order to 
obtain the transfer efficiency for the 4-level+ model at 
two-photon resonance, (8) = 0, we apply a Gaussian 
filter to this profile, centred at 8 — and then integrate 
over the result. By performing this integration process 
for a variety of delays we obtain the 4-levcl+ delay scan 
plotted here in Fig. [5J 

In the situation in Fig. [9] the Raman inhomogeneous 
lincwidth is 50 kHz and the Rabi frequency on the strong 
transition is 1 MHz, with an associated weak Rabi fre- 
quency of 370 kHz. Thus for these large Rabi frequen- 
cies and small inhomogeneous broadening we expect to 
be in the previously well-studied situation where the in- 
homogeneous broadening can be ignored @, HH, as is 
demonstrated in these models. 

The oscillations seen in the 4-level models for posi- 
tive delays are thought to be due to the fact that the 
models were performed at A = 0, optical resonance. At 
optical resonance it has been shown that the population 
transferral via 6-STIRAP oscillates with the adiabacity 
parameter (28[, which in our case depends solely on the 
pulse delay between the pump and Stokes pulses. 

In Fig. [TU] is shown the one-photon linewidth for the 
various models along with that observed from experi- 
ment. We see in this figure that the 4-levcl+ model agrees 
closely with the experimentally obtained result for this 
particular ground state splitting and Rabi frequency. 



In Fig. [TT] we compare the experimentally obtained 
two-photon detuning scan with that of the different mod- 
els. In the conditions of Fig. [11] the width of the two- 
photon detuning profile is much larger than the inhomo- 
geneous broadening, of 77.9 kHz. As a consequence, the 
presence of the inhomogeneous line has little effect on the 
transfer efficiency. The averaged profile given by the 4- 
level+ model hardly departs from the plain 4-level model. 
The computed profiles closely agree with the experimen- 
tal data. In this experiment the strong Rabi frequency 
was 790 kHz, giving a Rabi frequency of 290 kHz on the 
weak transition, which is still much wider than the inho- 
mogeneous width. 
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FIG. 10: (colour online) The solid (black) line shows the 
experimentally determined one-photon lineshape. This line- 
shape is also determined with the standard 3-level model, dot- 
ted (purple) line, the standard 4-level model, dashed (red) 
line, and the 4-level model with the inhomogeneous linewidth 
of the hyperfine transition included, dash-dot (blue) line, la- 
belled as 4-level+. The Rabi frequency of the strong transi- 
tion was 510 kHz, the ground state splitting was 7.01 MHz 
and the inhomogeneous width of the Raman transition was 
77.9 kHz. 



FIG. 11: (colour online) The two-photon detuning value is 
varied and the transfer efficiency is recorded experimentally, 
black squares. The dotted (purple) line is the output from the 
standard 3-level model, the dashed (red) line is the output 
from the standard 4-level model and the dash-dot (blue) line 
is the output from the 4-level+ model, where in this case the 
two-photon detuning value plotted is actually the average two- 
photon detuning value. The ground state hyperfine splitting 
was 7.1 MHz with an inhomogeneous width of 77.9 kHz and 
a Rabi frequency on the strong transition of 788 kHz. 

We now follow with an investigation into the effect 
the inhomogeneous width of the Raman transition has 
on the maximum transfer efficiency and the two-photon 
lincwidth as we vary the applied Rabi frequency. 

Shown in Fig. [12] is the effect on the maximum transfer 
efficiency as the Rabi frequency is varied. It is at low Rabi 
frequencies that the inhomogeneously broadened Raman 
transition is expected to have the largest effect on the 
transfer efficiency and we see that this is matched with 
the experimental and numerical results. 

At low Rabi frequencies we expect a lower maximum 
transfer efficiency for the 4-level+ model when compared 
to that of the standard 4-level model. In the standard 
4-level model, at 8=0, all the ions are at two-photon res- 
onance. Thus the majority of the population transferral 
occurs via efficient adiabatic transitions. In the 4-lcvcl+ 
modcl, at (S)=0, only a small portion of the ions are 
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FIG. 12: (colour online) We vary the Rabi frequency whilst 
the inhomogeneous width and ground state splitting remain 
constant. The splitting is constant at 7.1 MHz and the 
width of the inhomogeneously broadened Raman transition is 
79 kHz. The Rabi frequencies shown are those of the strong 
transition in each case. The squares are the results obtained 
from experiments. The solid (black) line is the simulation 
from the standard 3-level model, the dashed (red) line is the 
result from the standard 4-level model whilst the dash-dot 
(blue) line is the 4-level+ model. 



at two-photon resonance. The rest will have a non-zero 
two-photon detuning value, dependant on their position 
in the inhomogeneous line. Thus the majority of popu- 
lation transfer in this model occurs via the less efficient 
diabatic transitions, reducing the maximum transfer effi- 
ciency. As the Rabi frequency is increased a greater pro- 
portion of the transitions are adiabatic in the 4-lcvcl+ 
modcl, thus steadily increasing the maximum transfer 
efficiency, until it matches that of the standard 4-level 
model. 

In this figure we also see a difference between the 3- 
level and 4-level models. It is accepted that in a three- 
level system the population transfer is complete for a 
sufficiently high Rabi frequency satisfying the adiabac- 
ity criteria of either Eqns. \5\ or [Jj The transfer process 
then simply becomes more robust as the Rabi frequency 
is increased. This is seen in Fig. [T^] where the transfer 
efficiency of the 3-level model rapidly reaches it's max- 
imum value of ~94% and then remains constant as the 
Rabi frequency is increased. The maximum transfer effi- 
ciency in this model is not 100% due to the inclusion of 
relaxation and decoherence decay terms in the model. 

The situation is different for both the 4-level models 
where the transfer efficiency slowly reaches it's maximum 
value before then decreasing for larger Rabi frequencies. 

This reduction in the transfer efficiency might be due 
to the "connectivity problem" , which has previously been 
observed in multi-level systems [2!| HO, HH . The connec- 
tivity problem means that whilst there may be an adiab- 
taic dressed state overlapping with the initial bare state 
at t — > — oo and another adiabatic dressed state overlap- 



ping with the desired final bare state at t — -> oo there is 
no adiabatic or diabatic pathway connecting these two 
dressed states, due to the additional energy level and 
the multiple excitations of the transitions, resulting in 
no population transferral. It has previously been demon- 
strated that for multi-level systems it is possible that the 
maximum transfer efficiency will reduce for large Rabi 
frequencies 0, [3(| due to this problem. It is also possi- 
ble that at large Rabi frequencies the Stokes pulse could 
be strongly interacting with populated levels, resulting 
in optical pumping rather than STIRAP processes and 
thus reducing the transfer efficiency. 

Both experimental results and simulations show that 
the transfer efficiency, below that of maximum efficiency, 
is highly dependant on the Rabi frequency. From the adi- 
abacity condition of Eqn.[7]we see that for a constant one- 
photon detuning and pulse duration that it is the varia- 
tion of the Rabi frequency which will affect the transfer 
efficiency. As the Rabi frequency is reduced a greater 
proportion of the population transferral will occur via 
non-adiabatic transitions, thus reducing the transfer effi- 
ciency due to the associated radiative losses. 

We see explicitly in Fig.[T2]the dependence of the width 
of the two-photon detuning scan on the Rabi frequency. 
We see that as we reduce the Rabi frequency, both the ex- 
perimentally obtained maximum transfer efficiency and 
the two-photon detuning width is decreased. 









strong 


* 


950kHz 


♦ 


790kHz 




500kHz 


• 


370kHz 


■ 


200kHz 




-200 200 400 

Two-Photon Detuning (kHz) 



FIG. 13: (colour online) We varied the Rabi frequency for 
a constant ground state splitting whilst scanning the two- 
photon detuning. The symbols are the experimentally ob- 
tained results whilst the lines are the fitted Gaussians. The 
Rabi frequencies of the strong transition are listed for each 
scan. The ground state splitting in each case was 7.3 MHz 
with an associated inhomogeneous width of the Raman tran- 
sition of 80 kHz. 

Reducing the effective Rabi frequency will affect the 
maximum transfer efficiency in two ways. In order to 
achieve the maximum possible transfer efficiency the 
weak Rabi frequency should be much larger than the in- 
homogeneous width of the Raman transition and at the 
same time the adiabacity condition of Eqn. [7] must be 
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satisfied. 

Away from two-photon resonance population transfer- 
ral is only possible via diabatic transitions. As the Rabi 
frequency is reduced these diabatic transitions become 
weaker and hence it is thought that the weakening of 
these transitions reduces the two-photon width with de- 
creasing Rabi frequency. 

We monitored in Fig. [H] the two-photon detuning 
width as we varied both the Rabi frequency and the in- 
homogeneous width of the system. We found that the 
two-photon width was highly dependent on the Rabi fre- 
quency. It was also dependent on the inhomogeneous 
width of the hyperfine transition, but less so and only at 
low Rabi frequencies. 
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at low Rabi frequencies is to be expected. When the av- 
erage two-photon detuning value is varied slightly in the 
4-level+ model there can still be some ions at two-photon 
resonance, though less than if the average two-photon de- 
tuning was zero. Whilst in the standard 4-level model, 
once there is a non-zero two-photon detuning value there 
are no ions at two-photon resonance. Therefore there is 
greater population transfer in the 4-lcvcl+ model for non- 
zero two-photon detuning values, and hence a broader 
two-photon width. This effect is larger for the model 
with larger Raman inhomogeneous broadening as is to 
be expected. 

As the Rabi frequency is increased the two-photon 
linewidths of the various models converge. It is seen 
that the 4-lcvel+ model with the smaller inhomogeneous 
width, r = 80 kHz, converges at a smaller Rabi fre- 
quency than the model with a larger inhomogeneous 
width, r = 210 kHz, where T represnets the inhomo- 
geneous linewidth. We see that in each case the models 
with inhomogeneous broadening converge for a Rabi fre- 
quency of the strong transition of approximately six times 
the inhomogeneous width. 

This effect can also be seen in the experimental val- 
ues where we see that for low Rabi frequencies the two- 
photon detuning widths appear to be dependent on the 
Raman inhomogeneous widths. This effect appears to 
decrease as the Rabi frequency is increased. 



FIG. 14: (colour online) We varied the applied Rabi frequency 
and recorded the width of the two-photon detuning scan for 
different applied magnetic fields, and hence different Raman 
inhomogeneous widths. The associated Raman inhomoge- 
neous widths for the different ground state splitting values 
are F ~80 kHz for A 9 =7 MHz, F ~140 kHz for A 9 =13 MHz 
and T ~210 kHz for A 9 =19 MHz, where T is the width of the 
inhomogeneously broadened Raman transition. The exper- 
imental two-photon detuning values, shown as the symbols, 
are the FWHM of a fitted Gaussian to the experimental scans. 
The dashed (red) line is the standard 4-level model linewidths 
for a ground state splitting of A 9 =7 MHz. The dash-dot 
(blue) and the solid (black) lines are the linewidths from the 
4-level+ model for ground state splittings of A 9 =7 MHz and 
19 MHz respectively. 

It was found numerically that there was a negligible 
difference between the two-photon linewidths from the 
standard 4-level models for different ground state split- 
tings, thus only the result from a splitting of A g =7 MHz 
is displayed here. Once the inhomogeneous width is in- 
cluded in the model we see that the two-photon width 
is dependent on the inhomogeneous width at low Rabi 
frequencies. It is seen in this figure that the larger the 
inhomogeneous broadening and the smaller the Rabi fre- 
quency, the broader the two-photon linewidth of the 4- 
level+ model is when compared to that of the standard 
4-level model. 

The fact that the two-photon width is broader for the 
4-level+ model over the standard 4-level model, at least 



VI. CONCLUSIONS 

We have investigated STIRAP in Tm 3+ : YAG both ex- 
perimentally and numerically. The experimental demon- 
strations presented here were the first such demonstra- 
tions of STIRAP in Tm 3 +:YAG. We observed robust 
population transfer with a maximum transfer efficiency 
of -90% for STIRAP. We were also able to observe b- 
STIRAP in this system due to the long population life- 
time of the excited level as compared to the interaction 
time of the experiment. Maximum transfer efficiencies of 
—45% were observed for 6-STIRAP, significantly larger 
than the transfer efficiency for coincident Stokes and 
pump pulses. 

We were also able to completely model the Tm 3+ :YAG 
system as this simple system consists of only 4 energy lev- 
els and a metastable state in a magnetic field. We found 
that it was necessary to include the second excited state 
hyperfine level and the multiple excitations of the optical 
transitions in order for our numerical model to match the 
experimental results. It was also seen that the inclusion 
of the Raman inhomogeneous broadening to our standard 
4-level model altered the output and gave us simulated 
results more closely matched with those obtained exper- 
imentally, particularly at low Rabi frequencies. 

We also discovered that the two-photon linewidth is 
broader when the inhomogeneous broadening is intro- 
duced. This was thought to be due to more efficient pop- 
ulation transferral at non-zero two-photon detuning val- 
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ues with the inhomogcncous broadening present. It was 
seen that the two-photon linewidth determined from the 
model with inhomogcncous broadening converged with 
that of the standard 4-level model as the Rabi frequency 
was increased and the effect of the inhomogeneous broad- 
ening was reduced. 



It was also found that at low Rabi frequencies the 
presence of Raman inhomogeneous broadening reduced 
the maximum obtainable transfer efficiency. But at high 
Rabi frequencies, ie. frequencies ~10 times the inhomo- 
geneous width, the inhomogeneous width made no differ- 
ence to the obtained maximum transfer efficiency. 
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